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Effectiveness of a second phase in metallic glass heterostructures to improve mechanical 
properties varies widely. Unfortunately, methods to fabricate such heterostructures like foams 
and composites do not allow controlled variation of structural features. Here we report a 
novel strategy, which allows us to vary heterostructural features independently, thereby 
enabling a systematic and quantitative study. Our approach reveals the optimal micro- 
structural architecture for metallic glass heterostructures to achieve tensile ductility. Critical 
design aspect is a soft second phase, which is most effective when spacing between the 
second phase assumes the critical crack length of the metallic glass. This spacing should 
coincide with the second phase's size, and beyond, the specific second phase morphology of 
the heterostructure is crucial. These toughening strategies are only effective in samples that 
are large compared with the spacing of the second phase. The identified design aspects 
provide guidance in designing tensile ductility into metallic glasses. 
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Despite high strength and elasticity combined with plastic- 
like processing ability, wide spread proliferation of 
metallic glasses (MGs) in structural applications has been 
stymied by their lack of tensile ductility 1 . Upon yielding, shear 
strain in MGs is highly localized into narrow, ~ 10-nm-wide 
bands. The development of such shear bands in MGs does not 
necessarily result in fracture. For example, when shear bands are 
spatially confined under certain loads and in small sizes, global 
plasticity enabled by the formation of a large number of shear 
bands has been observed 2-15 . In compression and bending, 
multiple shear band formation without crack formation has been 
reported 5 ' 6 ' 11 ' 13 ' 16-18 . In nano-sized samples, it has been argued 
through a Griffith-like criterion that the elastic energy release 
does no longer offset the energy barrier to form a shear band 10 . 
Hence, on such a small scale, a deformation mechanism not 
relying on shear bands can result in tensile ductility 8 ' 10 ' 12 ' 19 . 

One emerging strategy to couple the attractive properties of 
MGs with plasticity is to introduce a second phase, which reflects 
and/or absorbs shear bands 2 ' 4 ' 6 ' 9 ' 20-25 . It has been suggested 
that this strategy should be particularly successful when spacing 
of the second phase coincides with the critical crack length of the 
MG matrix phase 9 . However, verifying this hypothesis as the 
origin of the toughening mechanism has been challenging as 
fabrication methods of MG composites and foams are complex, 
and do not allow independent and systematic variation of 
microstructural features, such as phase spacing, size, shape and 
volume fraction. 

Here, we report a novel approach that allows us to completely 
independently vary microstructural features with high precision, 
and determine their effect on mechanical properties. Using this 
systematic approach, we precisely determine that ductility 
increases with increasing the soft second phase spacing up to a 
spacing comparable to the critical crack length of the matrix 
phase. Furthermore, we emphasize the importance of morpho- 
logy, and identify and quantify the ideal heterostructure, which 
maximizes strength and toughness. The fracture strength of the 
optimized MG heterostructure almost reaches to the fracture 
strength of the monolithic MG while fracture strain and overall 
toughness of the MG heterostructure become 3.75 and 5 times 



higher than that of the monolithic MG, respectively. Our 
approach reveals the microscopic mechanism that results in 
tensile ductility of MG heterostructures, and provides guidance to 
design tensile ductility into MGs. 

Results 

Fabrication of MG heterostructures. Our method is composed 
of a two-step fabrication method, fabrication of a mould and 
replication of this mould with a MG (Fig. la). Silicon lithography 
is used to fabricate the moulds, and subsequent deep reactive ion 
etching allows to achieve deep cavities in the mould in excess of 
400 Jim for the required sample depth. Fabricated moulds com- 
prise the (negative) heterostructure, which is embedded into the 
gauge section of the tensile test samples. The moulds are repli- 
cated with Zr 3 5Ti3 0 Cu7.5Be 2 7.5 MG former through thermoplastic 
forming-based compression moulding 26-29 , and the fully 
amorphous state of these heterostructures is subsequently 
confirmed by differential scanning calorimetry (DSC) and X-ray 
diffraction (XRD). Noteworthy is the precision of this technique 
(Fig. lb), where features of the CAD drawing are replicated into 
the MG heterostructure through the Si mould with a precision 
better than 1 jam. 



Mechanical characterization of MG heterostructures. Tensile 
stress-strain curves for various MG heterostructures, and as a 
comparison, for a monolithic MG sample, are shown in Fig. 2a. 
The monolithic MG fabricated with our method exhibits an 
elastic strain limit of 2% and yield strength of 1,750 MPa, which is 
almost identical to the elastic strain limit of 2% and yield strength 
of 1,765 MPa, typically measured in the as-cast state of bulk 
monolithic rod-shaped samples in tension 30 . These almost 
identical values confirm that our fabrication method (Fig. 1) 
itself does not affect the mechanical properties. MG hetero- 
structures with specific microstructural architecture exhibit a 
remarkably different mechanical response. For a microstructural 
architecture with pore diameter (d) and spacing (s) of 50 urn, the 
fracture strain increases dramatically to &^7.5%. At the same 
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Figure 1 | Fabrication steps of MG heterostructures. (a) MG (in other words, bulk metallic glasses -BMGs) heterostructures are synthesized through are 
synthesized through Si lithography to create negative moulds. Such moulds are subsequently replicated through thermoplastic forming (TPF)-based 
compression moulding. Excessive material is removed by sanding/polishing/thermoplastic hot scraping, and the sample is released from the mould by 
selective KOH etching, (b) MG heterostructure realized in a tensile test sample, and close-up image of its gauge section (scale bar 3 mm). Comparison of 
the Si mould with the original CAD drawing reveals the control and precision of our process (scale bars 100 |im). 
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time, Young's Modulus of the heterostructure decreases, whereas 
the nominal fracture strength, cr f , remains relatively unaltered. 

In order to deconvulate elastic and plastic deformation, we 
carried out cyclic loading with gradually increasing loads 
(Fig. 2b). This figure shows the contribution of elastic/plastic 
components of MG heterostructures until fracture under uniaxial 
tension, which changes with second phase morphology. A 
heterostructure with d=\00 Jim and s = 50 Jim was selected over 
d — s — 50[im heterostructure for sample preparation conveni- 
ence. We found that of the 5.5% strain to failure, -70% (3.9% 
overall) corresponds to elastic and ~30% 
(1.6% overall) corresponds to plastic deformation. The large 
elastic deformation originates from the AB pore- stacking 
order. According to cyclic results with different d and s (not 
shown here), the elastic to plastic deformation ratio remains 
relatively unchanged, which is estimated to show around 2.6% 
plastic and 4.8% elastic deformation for d — s — 50 um hetero- 
structure. A heterostructure with the same pore diameter and 
spacing (d=100um, s = 50jim) but AA stacking exhibits a 
significantly lower strain to failure of 2.6%, where only ~40% 
(1.1% overall) deformation is elastic (Fig. 2c). Therefore, the 
geometrical effect caused by AB stacking increases the elasticity, 
and as a result, enhances the plastic deformation of the MG 
heterostructure. 

When comparing MG heterostructures with identical pore size, 
but increasing pore spacing from 25-100 Jim, both cr f and Sf 
increases. However, when the spacing exceeds 100 Jim, both 
values drop dramatically (Fig. 2d). The highest performance in 
terms of <7 f and e f is present in heterostructures with a spacing 
between 100 Jim and 200 urn. This spacing is close to the critical 
crack length, a c & 120-260 um, of Zr 3 5Ti3 0 Cu7. 5 Be27.5, which is 
calculated according to a (1/ti) (K c /<j) 2 (dimensionless factor is 
taken as 1 for circular pores), with a — 1,750 MPa, and K c ~35 
MPam 1/2 (unpublished data)-50 MPam 1/2 (ref. 31). Precise data 



for K c are rare for MGs, which is believed to originate from 
difficulties in sample preparation and a high notch radius 
sensitivity 32 . 



Deformation behaviour as a function of d/s. In order to 
understand the microscopic origin of the MG heterostructures' 
tensile ductility, we carried out scanning electron microscopy 
imaging (Fig. 3). For the d — s = 50jim heterostructure, which 
was deformed until failure (e f ^7.5%), a large number of shear 
bands form between the pores throughout the entire gauge sec- 
tion (Fig. 3a). These shear bands carry the plastic deformation, 
and their large number results in observed global tensile ductility. 
Shear bands initiate at pores, where they follow the direction 
towards the neighbouring pore, in a similar pattern as observed in 
double notched MG samples 33 (Fig. 3b). Finite element modelling 
(FEM) simulations show that, as deformation proceeds, the 
complex stress field (Fig. 3c) redirects shear bands towards 
neighbouring pores. Fracture occurs typically along the shortest 
distance between pores, which is for the considered sample 
dimensions, perpendicular to the uniaxial loading direction 
(Fig. 3d). At the edges of the sample, fracture chooses a 
different path by following the shortest distance in the diagonal 
pore direction (inset in Fig. 3d). 

It should be mentioned that the magnitude of the stress 
concentrations depends also on the shear modulus of the second 
phase, where a second phase with significantly different shear 
modulus than the matrix causes high stress concentrations. From 
the FEM results, we conclude that the lowest stress concentrations 
are present for a second phase with a modulus similar to the 
matrix. However, in this situation, elastic strain values as shown 
in Fig. 2b cannot be achieved. 

The experimental results (Fig. 2a,d) reveal the pores' diameter 
to spacing ratio (d/s) as the critical factor that controls the 
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Figure 2 | Quasi-static tensile characterization of MG heterostructures. Nominal fracture strength of the samples is calculated from the effective area, 
(a) Heterostructures with various pore sizes compared with monolithic MG. Increased cyclic loading of MG heterostructures with d/s = 2. (b) AB stacking 
reveals 3.9% elastic deformation and 1.6% plastic deformation, (c) The amount of deformation in AA stacking is 1.1% elastic and 1.5% plastic, (d) Tensile 
ductility and strength increases with increasing pore spacing until this spacing exceeds the plastic zone size. 
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Figure 3 | Microscopic analysis of the deformation mechanism in MG heterostructures. (a) Global tensile ductility originates from the formation of 
multiple shear bands throughout the sample (scale bar, 50 (im). (b,c) Shear concentrations at the pores initiate the formation of shear bands, and 
redirect their propagation towards the nearest neighbour pore (scale bar, 100 |im). (d) Fracture of the MG heterostructure occurs perpendicular to the 
loading direction (scale bar, 1 mm). Fracture is redirected at the sample's edges, following the shortest distance between the pores (inset— scale bar, 
100 urn). 



mechanical behaviour of MG heterostructures (Fig. 4a). Fracture 
strength decreases with increasing d/s. Fracture strain, on the 
other hand, does not increase monotonously, but instead, exhibits 
a maximum at d/s « 1 . The maximum in e f originates from an 
increase in the elasticity and plasticity of the AB- stacking 
heterostructure with increasing d/s, which is eventually over- 
compensated by rapidly increasing stress concentration between 
pores. As a consequence, the toughness of the MG hetero- 
structure ( (7t = Jq od&) exhibits a maximum at d/s ^ 1 of 
[7 T ^80MJm -3 , which is ~5 times higher than for the 
monolithic Zr 3 5Ti3 0 Cu7. 5 Be27.5 MG of l/ T =16MJm -3 
(Fig. 4b). It is noteworthy to state that the change in d and s 
have negligible effect on fracture stress and strain for the same 
d/s — 1 sample. Therefore, we concluded that d/s is the prominent 
factor in controlling the stress distribution between the pores, and 
thereby, the fracture stress and strain. 

To reveal the microscopic origin of the dependence of U T , Sf 
and <7f on d/s, we investigated the shear band propagation and 
distribution across the pores. Specifically, we characterized shear 
band patterns for different d/s ratio, and compared those with 
stress fields, which we determined through FEM. For d/s= 1, a 
large number of shear bands form throughout the sample 
(Fig. 4c). FEM analysis indicates that resultant von Mises stress in 
the middle region between the pores perpendicular to the loading 
direction of ~ 800 MPa is lower than the applied stress of 
1,000 MPa. Stress concentrations are confined to small, non- 
percolating regions adjacent to the pores (Fig. 4d). For d/s = 4, 
only a very small number of shear bands form before fracture 
(Fig. 4f). Thus, the resulting toughness of the heterostructure is 
low (Fig. 4b). FEM results reveal that resultant von Mises stresses 
for d/s — 4 are significantly higher ( ~ 3,200 MPa) for an applied 



stress of 1,000 MPa, and stress concentrations overlap between 
pores (Fig. 4g). As a consequence, MG heterostructures with 
d/s — 4 yield at a lower applied load level, and hence, the total 
amount of plastic deformation until fracture is much smaller 
(Fig. 2d). As a result, only a small number of shear bands form 
before they develop into cracks. 



Effect of embrittlement and crystallization. In addition studying 
the effect of pore spacing on the MG-heterostructures' perfor- 
mance, we also varied the properties of the matrix material, 
and thereby, its critical crack length. An effective strategy to 
vary the critical crack length is to structurally relax (embrittle) the 
MG through thermal annealing 1 1,13,34_39 . We fabricated tensile 
samples out of Zr 35 Ti3 0 Cu 7 5 Be 2 7.5 MG former with d/s— 1.5, 
which were annealed at 558 K, ~25K below its calorimetric 
glass transition temperature. The degree of structural relaxation 
is determined through enthalpy recovery 40 (not shown here), 
and its effect on mechanical behaviour through bending tests 
(Fig. 5b). Tensile characterization shows that strain to failure 
and U T are dramatically affected by annealing (decrease in a c ) 
and crystallization (Fig. 5a,b). The completely structurally relaxed 
MGs for 36 h and longer, as well as the crystallized 
MG former approximate ideal brittle behaviour with 
K IC <5MPa m 1 ^ (refs 41,42), and exhibit a critical crack 
length below 1 um. In such MG heterostructures, a c < s, and as 
a consequence, the overall strain decreases to the elastic limit of 
~ 2%. Even though K IC is not directly quantified in this work, it 
can be concluded from Fig. 5b that for the as -cast and 28 h 
annealed heterostructures a c >s, hence shows large tensile 
ductility. 
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Figure 4 | Effect and origin of d/s on mechanical properties in MG heterostructures. (a) Fracture stress and fracture strain as a function of the ratio of 
diameter to spacing of the pores (d/s). (b) Toughness (L/ T ) versus d/s exhibits a pronounced maximum at d/s^1. Error bars represent the s.d. of the mean 
values from at least three test samples, (c) Corresponding scanning electron microscopy image of MG heterostructures with d/s = 1 showing multiple 
shear bands (scale bar, 50 (im). (d) For an applied average stress of 1,000 MPa, the spatially confined stress concentrations between the pores 
perpendicular to the loading direction are lower compared with heterostructures with d/s = 4, hence multiple shear bands form without crack formation, 
(e-g) Shear bands initiate at the shear stress concentration regions adjacent to the pores, (f-h) For d/s = 4, stress concentration pathways between the 
pores already exceed cr y for the same average applied stress, resulting in early fracture from only a few shear bands (scale bar, 50 |im). 
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Figure 5 | Influence of sub-T g annealing and crystallization on the ductility of heterostructures. (a) For samples that are completely structurally relaxed 
at 558 K, as well as for crystallized samples, ductility decreased dramatically, (b) Bending test reveals a significant decrease in plasticity for annealing times 
exceeding 28 h at 558 K. 



Discussion 

The aim of our approach is to identify and understand the 
mechanisms that result in tensile ductility in MG heterostructures 
such as foams, cellular heterostructures and, to a limited extent, 
composites. The strength and purpose of our method is the ability 



to precisely vary individual microstructural features completely 
independent, and determine the effect of such isolated variations 
on the mechanical response. As such, our method is not intended 
as a scalable method to fabricate MG heterostructures for 
applications, but as a method to identify and understand the 
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effectiveness of some MG heterostructures to guide design and 
synthesis of MG foams, cellular structures and composites. 

Our results revealed three critical aspects of a MG hetero- 
structure, which control mechanical behaviour. One is the 
difference between the second phase spacing and the critical 
crack length of the MG, s — a c . For s — a c < 0, shear bands 
that form between pores do not grow to a length where they 
transform into a crack. As a consequence, global tensile ductility 
can be achieved through the formation of multiple shear 
bands. For s — a c > 0, a shear band can develop into a crack 
before strain energy releases into the pore and, consequently, only 
one or a few shear bands form before crack formation, which 
impairs tensile ductility. Our finding has previously been broadly 
confirmed in MG heterostructures with a soft second phase, that 
is, G 2n d < G matrix , such as foams, cellular structures and 
composites 2 ' 9,20 ' 21,23-25 ; nevertheless, up to date, it has not been 
proven owing to experimental limitations. 

However, the s — a c criterion alone is insufficient in predicting 
the performance of MG heterostructures with a soft second phase 
owing to the effect of stress concentrations on the deformation 
mechanism. In addition to s — a c , we identified the ratio of second 
phase size to spacing as a critical feature for toughness and tensile 
ductility optimization. Maximum toughness is achieved for 
d/sttl, which corresponds to a density of ~50%. This value for 
d/s exhibits the best compromise between enhanced stress 
concentrations, which limits plastic deformation, and thin 
interporous regions, which enhance plasticity in cellular structures. 

Our design criteria is suggested to be readily applicable to 
cellular MGs, whereas the applicability for the effective design of 
MG composites is limited, which is owing to the complex nature 
of MG composites, and the mechanisms that toughen them. 
Experimental findings suggest influences on the toughening 
mechanism of shape and volume fraction of the second phase, 
interfacial strength and toughness, softness of the second phase 
reflected in shear modulus, and even dislocation and twinning in 
the second phase. Our finding of the drastically different 
mechanical response between MG heterostructures with AB and 
AA stacking exemplifies the importance of the second phase 
morphology and the overall complex nature of the toughening 
mechanism in MG composites. 

Furthermore, the softness of the second phase has been broadly 
suggested to have a crucial role in the toughening mechanism . 
Even though the pores in the MG heterostructure have a smaller 
shear modulus than the matrix, which has been recognized to be a 
common feature of effective toughening mechanisms, it is not 
clear if pores can be considered as a soft second phase, where 
their shear modulus is significantly lower than G matrix , and they 
lack an interface. 

We argue that the effectiveness of MG heterostructures also 
depends on the size (cross-section width perpendicular to the 
loading direction) of the sample, w, in comparison with s. Such a 
size effect originates from the fact that shear bands do not 
support tensile forces during the actual shearing event. As the 
cross-sectional area supporting a tensile load in samples, where w 
is not much larger than s, is significantly reduced when one shear 
band forms, the applied load results in exceeding a y throughout 
the sample. Therefore, small samples with w comparable to s fail 
catastrophically along one shear band. In large samples, w»s 
(assuming s — a c <0), upon formation of a shear band between 
two neighbouring pores, stresses between most of the other pores 
typically do not exceed cr y as the stress increases with the area 
reduction, w-s, owing to the initial shearing event being small. 
Our finding also explains why all reported MG composites, which 
are exhibiting tensile ductility, obey w » s. 

Our results also reveal the possibility to design elasticity 
through the elasticity and the modulus of the second phase in MG 
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heterostructures. As shown in Fig. 2b, we have demonstrated up 
to 3.9% of elastic strain for the heterostructure with d— 100 Jim 
and d/s — 2, and an AB- stacking order. Morphology in MG foams 
and composites is AB-like; however, in MG foams and 
composites, an increase in elasticity compared with the mono- 
lithic MG has not been reported yet. For MG composites, this is 
owing to the properties of the second phase, which exhibits low 
elasticity (typically below 0.5% and comparable modulus to the 
MG matrix), thereby preventing increased heterostructure 
elasticity. MG foams with controlled features have not been 
tested in tension yet. 

In summary, we identify the most effective heterostructure to 
result in tensile ductility and toughness in MGs. Such structure 
features a spacing of the second phase, which coincides with the 
critical crack length of the MG and with the second phase size. In 
addition, we find that heterostructure morphology is crucial and 
the size of the sample should be large compared with the second 
phases' spacing. 

Methods 

MG alloy production. Zr 3 5Ti3oCu7. 5 Be27.5 alloy was prepared by melting the ele- 
ments with purity better than 99.99% in an arc melter. To prepare an amorphous 
rod, the ingot was sealed in a quartz tube under vacuum, which was then heated 
above its liquidus temperature for 5 min, and subsequently quenched in water. 
X-ray diffractometer (XRD-6000 Shimadzu) and differential scanning calorimeter 
(Perkin Elmer Diamond DSC) confirmed the amorphous nature of the Zr-MG 
rods. 

Si mould fabrication. We design MG heterostructure patterns using AutoCAD 
201 1 and Layout Editor. The drawing layout was transferred to a photomask maker 
(Heidelberg DWL-66 Laser Mask Writer), and the photomask is generated using 
direct laser beam writing. A thin layer of photoresist (8-um AZ 5214 coating) on 
650-um-thick Si wafer was patterned by conventional photolithography using EVG 
620 Contact/Proximity Mask Aligner. The etch cycle for DRIE processing (Bosch 
etch) was carried out under 12 s, 130-s.c.c.m. SF 6 , 13-s.c.c.m. 0 2 , 600-W coil power, 
12-W platen power and 26 mtorr pressure, and the passivation cycle was performed 
under 8 s, 85-s.c.c.m. C 4 F 8 , 600-W coil power, 0-W platen power and 15 mtorr 
pressure. The overall etch rate was 4ummin - l , where the final depth of the Si 
moulds was measured to be around 400 ± 20 um. Residual resist was subsequently 
stripped from the surface using an etchant (PRS 3000), and cleaned by soaking first 
in acetone, then in isopropanol and finally in deionized water. The Si wafer was 
finally dried with N 2 under a pressure of 2 x 10 5 Pa. 

MG heterostructure fabrication. MG discs cut from the amorphous cast rods 
were pre-pressed to a thickness of around 800 um using custom-built heating 
system in conjunction with Instron 5569 tensile testing machine (50 kN maximum 
load capacity). Using thermoplastic-based compression moulding, MG hetero- 
structures were fabricated from Zr 3 5Ti3 0 Cu7. 5 Be27.5 MG former at 425° C for 60 s 
under a pressure of 50 MPa using the above described setup. The extra MG layer 
remaining on the top surface of the silicon was later removed by grinding and 
polishing using Buehler Metaserv 250 Grinder- Polisher. The MG tensile hetero- 
structures are subsequently taken out of the mould by Si mould etching in KOH for 
60 min at 100 °C. The amorphous heterostructure of the MG honeycomb was later 
confirmed by thermal and structural analysis using DSC and XRD, respectively. 
Mechanical characterization of the MG heterostructures was conducted using 
Instron 5543 Tensile Tester (lkN maximum load capacity) under quasi-static 
conditions (strain rate of 0.001 s _1 ). Simultaneous image recording is conducted 
by Keyence Digital Microscope VHX-500F. Scanning electron microscopy imaging 
(Model: Hitachi SU-70) was subsequently utilized to determine the fracture pattern 
and the distribution of shear bands between pores. 

Finite element simulation. COMSOL Multiphysics Engineering Simulation 
Software was utilized to conduct structural element analysis for AB- stacking MG 
heterostructures. Linear elastic material model was used to approximate the 
mechanical behaviour of the Zr-MGs. The thickness of the specimen was set to be 
300 um. The surface was divided into refined triangular meshes, and an applied 
stress of 1,000 MPa is evenly distributed among the pores. Von Mises stress, shear 
stress and elastic strain around the pores were measured and compared with AA- 
stacking MG heterostructure, and other AB-stacking MG heterostructures with a 
second phase such as polymer, soft metal and graphite. 

Bending test. Rectangular beams of 0.65 ± 0.05 mm thick were prepared from the 
standard samples. All sides of the beams were mirror polished, and the corners 
were rounded to avoid stress concentration during bending tests. Beams were bent 
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over mandrels with decreasing radius from 30-5.5 mm with corresponding strains 
of 1.2-6% under constant compressive strain rate of 0.005 s _1 using Instron 
Testing Machine. This range spans the elastic to the plastic strain region, which 
enables to observe changes in the bending ductility with high accuracy. XRD 
spectrum and diffraction scanning calorimetry thermogram confirms structural 
differences between as-fabricated and sub- and above-T g annealed MG hetero- 
structures (not shown here). 



Heat treatment of MG heterostructures. Zr 35 Ti3 0 Cu7 . 5 Be 2 7. 5 MG hetero- 
structures were annealed for different time intervals at 285 °C for sub-Tg embrit- 
tlement, whereas the crystallization was achieved by heating the MG 
heterostructure to 420 °C (above-Tg) for 1 h using Vulcan Box Annealing Furnace 
(Model No: 3-550). Embrittlement can be studied in a wide range of MG alloys. 
DSC and XRD measurements were conducted to determine the effect of different 
annealing processes on thermal and structural properties. 
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